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ABSTRACT 

Wc develop an automated spectral synthesis technique for the estimation of mctallicities ([Fc/H]) and 
carbon abundances ([C/Fe]) for metal-poor stars, including Carbon- Enhanced Metal- Poor (CEMP) 
stars, for which other methods may prove insufficient. This technique, autoMOOG, is designed to 
operate on relatively strong features visible in even low- to medium-resolution spectra, yielding results 
comparable to much more telescope- intensive high-resolution studies. We validate this method by 
comparison with 913 stars which have existing high-resolution and low- to medium- resolution spectra, 
and that cover a wide range of stellar parameters. We find that at low mctallicities ([Fe/H] < —2.0), 
we successfully recover both the metallicity and carbon abundance, where possible, with an accuracy 
of ~ 0.20 dcx. At higher mctallicities, due to issues of continuum placement in spectral normalization 
done prior to the running of autoMOOG, a general underestimate of the overall metallicity of a star 
is seen, although the carbon abundance is still successfully recovered. As a result, this method is 
only recommended for use on samples of stars of known sufficiently low metallicity. For these low- 
metallicity stars, however, autoMOOG performs much more consistently and quickly than similar, 
existing techniques, which should allow for analyses of large samples of metal-poor stars in the near 
future. Steps to improve and correct the continuum placement difficulties are being pursued. 

Subject headings: Galaxy: halo - stars: abundances 
niques: spectroscopic 



stars: carbon - stars: Population III - tech- 



1. INTRODUCTION 



One of the most intriguing results fro m modern surveys 
of me tal-poor stars, such as the HK (IBeers et al.lll985l 
1992 ) and Hamburg/ESO surveys (HES: lChristlieb et alJ 



20011) . is the large fraction of Carbon-Enhanced Metal- 
Poor (CEMP; [C/Fc] > +1.0) stars identified amongst 
the lowest metallicity stars. The exact fraction of CEMP 
stars at low metallicity has been deba ted in the literature , 
and is variously reported a s ~10% dFrebel et alJl2006D. 
~15% (jCohen et al.l I2005D or -20% (jLucatello et al.l 
1200ft) of stars with [ Fe/H] < -2.0: ~20%-2 5% of stars 
with [Fe/H] < -2.5 (jMarsteller et al l I2005D ; or nearly 
40% of stars with [Fe/H] < -3.5 (|Beers fc Christliebl 
l2005h . The fraction of CEMP stars is expected, on 
theore tical grounds, to increase with declining metal- 
licity (jTumlinson I l2007f ). The evolutionary states of 
sample stars can also play a signific ant role, due to 
the mixing-induced dilution of carbon (jAoki et al.ll2006l : 
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iLucatello et al.ll2006fl . 

To better specify the changing CEMP fraction with 
metallicity, much larger samples of stars are required. 
Existing studies have generally been limited to (small) 
samples of low-metallicity stars for which high-resolution 
spectral a bundances have be en obtained (with the ex- 
ception of iFrebel et alJ l2006f h Such spectra arc time 
consuming to obtain and simply fail to exist in large 
enough samples. Lower-resolution spectra are available 
in sufficient quantity, but the techniques used to deter- 
mine abundances are more inconsistent across parameter 
space, and largely based on empiri cal calibrations (e.g., 
IBeers et alJl!999t iRossi et alJl2005D . 

Here we discuss the development of an automatic spec- 
tral synthesis fitting routine, autoMOOG, which enables 
rapid spectral fitting of the low- to medium-resolution 
spectra available in existing and future large samples 
of metal-poor stars. Papers in preparation will use this 
technique, and others, on a much larger sample of stars 
than has been studied to date, to determine the fre- 
quency of carbon enhancement as a function of metal- 
licity. autoMOOG will also be vital for identifying the 
most useful and interesting targets from large data sets 
for future high-resolution follow-up. 

This paper is outlined as follows. Section 2 presents a 
discussion of the current techniques and a description of 
the new method that has been developed. Error analysis, 
performed on the program using synthetic spectra, is pre- 
sented in Section 3. Section 4 considers validation of the 
autoMOOG approach using medium-resolution spectra 
of stars with available high-resolution determinations of 
[Fe/H] and [C/Fe]. The important question of the specifi- 
cation of upper limits on [C/Fe] for warmer stars (where 
the prominent C features are quite weak) is addressed 
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in Section 5. Section 6 presents brief conclusions and a 
discussion of future work anticipated for refinement of 
autoMOOG. 

2. METHODOLOGY 

For metal-poor stars with only low- to medium- 
resolution spectroscopy available, stellar parameters such 
as metallicity ([Fc/H]) and carbon abundance ([C/Fe]) 
have primarily been estimated from broadband colors 
and line indices measuring the strength of the Call K 
line an d the CH G band-the KP and GP i ndices, respec- 
tively (|Beers et al.lll999HRossi et al.ll2005l) . However, for 
CEMP stars, the presence of strong carbon features may 
affect the me asured color s of a star, in particular the 
B magnitude. iRossi et alj (|2005l ) have attempted to ac- 
count f or this by us i ng J — K colors instead of B — V 
used bv lBeers et al.1 (|1999f ). Strong carbon features may 
also affect the measurement of the continu um in a side 
band of the KP index (as first indicated by ICohen et al.l 
l2005h . Such effects could lead to an underestimation of 
the stellar metallicity. Very strong carbon features can 
also "overflow" into the sidebands used for the measure- 
ment of the GP index, confounding estimates of [C /Fe] . 

It is useful to explore different approaches that should 
be relatively immune to these effects. One such method 
involves the use of model stellar atmospheres to gener- 
ate synthetic spectra, so that detailed fits of a given ob- 
served spectrum may be made, rather than suffer the loss 
of information inherent in techniques based on line in- 
dices. We have explor ed using the spectrum analysis code 
MOOG (jSnedenlll973r i. a grid of ATLAS9 ID, LTE stellar 
atmosphere models and atomic and molecular line lists 
(inclu ding lines of CH, CN, and C2; ICastelli fc Kuruczl 
|2003[ ) to carry out this approach. We have developed a 
routine to automatically run MOOG and minimize the 
residuals of the resulting model spectra when compared 
against certain critical regions of the observed data. 

Using the li ne index metho ds in concert with other 
techniques fsee lLee et al.l2006l) , T e / / and log g values are 
adopted and initial estimates are determined for [Fe/H] 
and [C/Fc]. From our grid of models (with spacing of 100 
K in T e ff, 0.2 in log g, and 0.5 in [M/H], or overall solar- 
scaled metallicity), we select the model that most closely 
matches our starting parameters. This model is used 
within MOOG, where individual abundances are then 
varied to the nearest 0.05 dex to fit a full line profile. If 
the metallicity is varied significantly, the next appropri- 
ate model in the grid is used instead. Any element not 
currently being fit is given the model solar-scaled value. 
However, the user may specify other individual elements 
to vary as well, generally by a fixed offset from the scaled 
value. Thus, for example, it would be possible to use CH 
to determine [C/Fe], and then use that value to fit CN 
to determine [N/Fe]. 

As in the line index method, the Ca 11 K line at 3933 
A is used to solve for the metallicity of a star. At low- 
to medium-resolution, no iron lines are strong enough 
to be cleanly measured, so calcium is used as a proxy 
for metallicity. However, unlike the line index method, 
which empirically calibrates the strength of the KP index 
to the metallicity, here the Ca 11 K line is fit to obtain 
a calcium abundance. Since calcium is an a-element, 
an appropriate, metallicity-dependcnt offset, to account 
for possible a-enhancements, must be used when deter- 
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Fig. 1. — Example of the fits (line) produced by autoMOOG of 
the Ca II K region (left) and of the G band (right) compared to the 
data (triangles). The numbers in the right panels correspond to 
the T e ff, log g, [Fe/H], and [C/Fe], respectively, of the star. Note 
the higher resolution of the BM stars compared to the HERES 
sample, showing that a range of spectral resolutions can easily be 
fit. These stars also show the effect of temperature and abundances 
on the appearance of the Ca II K line and the G band, as well as 
the appearance of an upper limit determination. 

minin g the metallicity of a star. Following iBeers et al.l 
(fl999l) . we use an a-enhancement applicable for Milky 
Way stars of: 

f if [Fe/H] > 

[a/Fe] = <^ -0.267 * [Fe/H] if - 1.5 < [Fe/H] < (1) 
{ 0.4 if [Fe/H] < -1.5 

The initial guess for [C /Fe] is also included in our metal- 
licity fi t to address the effect pointed out bv lCohen et al.l 
(|2005] ). 

Using the new, more precisely determined metallicity, 
[C/Fe] can similarly be fit using the CH G band at 4304 
A. However, for warmer stars, those with a lower car- 
bon abundance, or those that are very metal-poor, the G 
band may be too weak to be detected. In these cases, an 
upper limit for [C/Fe] is determined instead. In theory, 
with an updated [C/Fe], it may be necessary to iterate 
on the fit to the Ca II K line. However, extensive testing 
revealed that the metallicity determined using the input 
[C/Fe] and that obtained from the refined [C/Fe] varied 
by negligible amounts (if at all), suggesting that such 
iteration is not required at this precision. 

Once this approach was worked out by fitting the abun- 
dances of a number of stars by hand, a program, auto- 
MOOG, was developed to iteratively vary abundances 
and automatically minimize the \ 2 of the synthetic spec- 
tra produced by MOOG. autoMOOG allows the user to 
select any wavelength region for the comparison. For the 
results presented here, a fit was determined using the 
ranges 3910-3960 A for the Ca n K line, and 4275-4330 
A for the CH G band. An example of the fits produced 
by autoMOOG can be seen in Figure Q] 

As part of the automation, autoMOOG also uses the 
region around the Ca II K line to estimate the observed 
radial velocity, which is used to align the data with the 
model spectra. The velocity precision is directly propor- 
tional to the spectral resolution, which for our test spec- 
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TABLE 1 

Abundance Errors, and Their Dependence on T e ff AND 
log g 



Bin 



[Fc/H] 
— i — 



[C/Fe] 
— r 



CTlnt "T bW °"log a Total °"Int CT T e // CT log9 °"Total 



Errors in Bins of T, 



eff 



4000 K 0.162 0.306 0.236 0.351 0.267 0.238 0.207 0.267 

5000 K 0.058 0.168 0.073 0.174 0.058 0.107 0.064 0.111 

6000 K 0.093 0.139 0.070 0.139 0.088 0.097 0.065 0.097 

7000 K 0.041 0.123 0.046 0.125 0.033 0.104 0.085 0.130 

Errors in Bins of log g 

O 0.124 0.172 0.112 0.172 0.153 0.153 0.119 0.153 

4.6 0.139 0.154 0.139 0.154 0.152 0.161 0.141 0.161 

5.0 0.180 0.254 0.172 0.254 0.174 0.168 0.143 0.174 

1 This error is the combined errors due to the intrinsic error 
and due to the uncertainty in the input effective temperature. 
Since the individual errors are non-Gaussian, there is no simple 
way to disentangle them. 2 This error is due to the combination 
of the intrinsic error and the uncertainty in the input surface 
gravity. 3 This is an estimate of the total error based on the 
individual errors. Since the individually measured errors are 
not separable, the total error cannot be calculated directly; this 
value is an approximation based on the measured values. 

tra (R ~ 4000) is typically <10 km s — 1 . Continuum sub- 
traction and spectral renormalization occur prior to the 
running of autoMOOG. Assuming this renormalization 
to be valid, autoMOOG applies a small additive offset to 
the local continuum level of the normalized, continuum- 
subtracted observed spectra such that it matches the lo- 
cal continuum level of the comparison synthetic spectra 
without altering the effective line strength. As we will 
see in Section 4, this assumption likely becomes less valid 
at higher metallicities, where the apparent continuum is 
significantly suppressed in the vicinity of the Ca II line, 
leading to incorrect continuum subtraction. 

3. ERROR ANALYSIS 

To estimate the errors on the derived abundances, a 
set of noiseless synthetic spectra across the grid of pa- 
rameters was generated. Since these spectra had known 
parameters, the abundance determinations could be di- 
rectly compared with the proper values. These spectra 
were first run through autoMOOG to determine any in- 
trinsic scatter in the estimates. These errors, and their 
rough dependence on T e ff and log g, are listed in Table 
[TJ The errors are roughly constant across both of these 
parameters, with the exception of very cool giants where 
the scatter increases. Also, note that the minimum er- 
rors are limited by the resolution of the abundance grid 
(0.05 dex). 

The intrinsic errors are nonzero due to the nature of the 
abundance determination. The program must first de- 
termine a cross-correlation velocity and continuum level 
in order to shift the data to match the model spectra. 
The finite resolution of the data and model spectra, as 
well as the finite grid of velocity shifts and input model 
abundances, occasionally leads to small offsets in abun- 
dance determinations. The uncertainty introduced by 
these various fitting steps, including the velocity and con- 
tinuum shift, as well as other sources of error, are then 
included in the observed intrinsic scatter. Note, however, 
that this scatter is still quite small. 

Since external values for T e f f and log g are adopted in 
choosing a model, appropriate errors in these parameters 
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Fig. 2.— Total error for [Fe/H] (solid line) and [C/Fe] (dotted 
line) vs. signal-to-noise ratio. 

(cr Te// = 200 K, criogg = 0.4: lLee et aT1l2008f ) need to be 
translated into errors on the abundances. Since [Fe/H] 
and [C/Fe] are fit internally, using the external values 
only as a starting point, the errors on these parameters 
do not need to be considered. 

To translate these errors, the same synthetic spectra 
used to determine the intrinsic error were refit with incor- 
rect stellar parameters (±100, ±200 K in T e ff, and ±0.2, 
±0.4 in log g, which correspond to the line index errors 
and are related to the spacing of the model grid) . The 
resulting combined errors for [Fe/H] and [C/Fe] based on 
the uncertainties in these parameters and the intrinsic er- 
rors are shown in Table[T] Total errors are then estimated 
from these individual sources by assuming approximate 
Gaussian distributions. Since the errors are essentially 
constant across log g, the total errors are most reason- 
ably read from the bins of varying temperature, and are 
thus less than ~ 0.20 dex for [Fe/H] and ~ 0.15 dex for 
[C/Fe], for stars warmer than 4000 K. 

The error in the abundances due to the changing signal 
to noise of the data also needs to be considered. To deter- 
mine this, varying amounts of noise were injected into the 
noise-free synthetic spectra, which were run through au- 
toMOOG again. This error combined with the previous 
errors (assuming crrpe/H] = 0-20 dex and o"[c/Fc] — 0.15 
dex) are shown in Figured] The total errors are then less 
than 0.25 dex for S/N > 10 for both [Fe/H] and [C/Fe]. 

4. VALIDATION 

Validation of the autoMOOG approach requires com- 
parison of the results, as obtained from medium- 
resolution spectra, with stars that have abundances de- 
termined from high- resolution analyses. Differences be- 
tween high- resolution and lower- resolution analyses, such 
as differences in adopted stellar parameters, different 
lines used for individual abundances, and differences in 
fitting techniques, may all lead to small scatter between 
results. However, lower-resolution techniques are often 
empirically calibrated to higher-resolution studies, which 
somewhat mitigates this scatter. Comparison of the 
autoMOOG technique with other methods designed for 
lower-resolution spectra, such as the line index method, 



4 



Marstcllcr ct al. 



- 

[x. -2 





i 


1 


i 

Bidelman & MacConnell □ 

HERES * ~ 


















h 

..* T* * = i nSe b" ** ' 


(a) 

1 h- 


1 1 — 


— 1 |— 


— 1 1 1 1 
























. 

* * * * t 




; -*i 




I 


i 





[Fe/H] h 



[Fe/H] . -1.75 
[Fe/H] . -0.60 




3920 3930 3940 3950 

Wavelength (A) 



Fig. 3. — Comparisons of the mctallicity determined from (a) 
autoMOOG and (b) the line index method to the high-resolution 
results. In both cases, the higher metallicity, carbon-normal BM 
stars are plotted as open squares, while the more metal-poor, oc- 
casionally carbon-rich HERES stars are plotted as filled triangles. 
The solid line corresponds to agreement with the various abun- 
dances. The dashed lines to either side correspond to a difference 
of ±0.5 dex. 

would also be valuable, particularly since they would be 
looking at the same features. These comparisons are per- 
formed for two sets of data, as described below. 

4.1. The Data 

The first set of stars analyzed are from the catalog 
of iBidelman fc MacConnelll (|1973l hereafter BM), orig- 
inally identified as metal-weak candidates in the south- 
ern sky. This sample ranges from around [Fe/H] = — 2.5 
up to solar mctallicity, and thus represents a modestly 
metal-poor, and predominantly carbon-normal, sample. 
The high-resolution values used here come from the liter- 
ature compilation of|C^reOiISilobere£^n ([2001). The 
medium-resolution (~ 1 A, S/N ~ 100) spectra that are 
being fit are described in T. C. Beers et al (2009), in 
preparation. 

The other data are a sample of metal-poor stars ob- 
served as part of the H amburg/ESO R-p rocess Enhanced 
Star survey (HERES: IChristlieb et aDl2004h . Initially 
selected from the full Hamburg/ESO and HK surveys, 
the HERES project obtained moderate-resolution (~ 2 
A) follow-up spectroscopy for metal-poor giant candi- 
dates. For bright stars confirmed to have [Fe/H] < —2.5, 
higher-resolution {R = 20,000, S/ N ~ 50, A = 3760 - 
4980 A) spectroscopy was obtained (|Barklem et al.|[2005t 
iLucatello etHI l2006t ) . Their goal was to identify the 
small fraction of these stars that showed large enhance- 
ments of r-process elements, so that investigations into 
the nature of this process could be conducted. As a con- 
sequence of these observations, several hundred metal- 
poor stars, many of them carbon-enhanced, were ob- 
served at moderate and high resolution. 

These then make ideal samples for the present anal- 
ysis, as the medium-resolution data can be run through 
autoMOOG, and the high-resolution results are available 
for the purposes of comparison. 



Fig. 4. — Synthetic spectrum near the Ca II K line assuming 
[Fe/H] = -1.75 (dashed line) and [Fe/H] = -0.60 (solid line) for 
a star of T e ff = 4800, log g = 2.2. Overplotted are points cor- 
responding to data from a star in the BM sample, using the data 
created after the continuum fit (filled triangles) and that same data 
scaled to match the solid line (open triangles). Note the different 
apparent continuum levels for the metal-poor spectrum vs. the 
metal-rich spectrum. 

4.2. Results 

Due to the high quality spectra available for the major- 
ity of the medium-resolution data, metallicities and car- 
bon abundances were determined for nearly all the stars 
in the samples. A total of 913 stars were fit, including 
529 from the BM sample and an additional 384 from the 
HERES sample. In some cases, stars had no metallicity 
determination available from the line index method. In 
addition, many of the BM stars lacked high-resolution 
metallicities, and almost none had high-resolution car- 
bon abundances. In all of these cases, abundances were 
obtained using autoMOOG, but obviously no comparison 
was possible. 

When comparing autoMOOG results for metallicity 
with the high-resolution analysis (Figure^), it can be 
seen that there exists a generally good agreement at low 
metallicities, with a small scatter and no appreciable off- 
set. Due to issues of continuum placement in the normal- 
ization of higher metallicity ([Fe/H] > —2.0) input spec- 
tra, an increasing underestimate of metallicity is seen. 

During the process of continuum subtraction and spec- 
tral normalization prior to the running of autoMOOG, 
a continuum is placed along the apparent continuum of 
a stellar spectrum. For high-metallicity stars observed 
at low spectral resolution, the multitude of metallic lines 
causes a suppression in the apparent continuum, mask- 
ing the true continuum level. Continuum subtraction 
techniques will improperly normalize such spectra, and 
subsequent lines fit within autoMOOG, such as calcium, 
will then appear weaker than they truly are, resulting in 
a general underestimate in the metallicity estimates. 

This effect is demonstrated in Figured] Synthetic spec- 
tra are plotted for models with T e ff = 4800 K and log 
g = 2.2, with possible metallicities of [Fe/H] = —1.75 
(dashed line) and [Fc/H] = —0.60 (solid line). Overplot- 
ted on these synthetic spectra are the data from star 
CD-57:949, taken from the BM sample. The filled tri- 
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angles correspond to the flattened, normalized spectrum 
as produced from the continuum fit to the data. The 
open triangles, on the other hand, are the same data af- 
ter application of a scaling factor selected to best match 
the higher metallicity synthetic spectra. Based on the 
input, filled triangle spectra, autoMOOG appropriately 
selected a metallicity of [Fc/H] = —1.75 for this star. 
The fit is far from perfect, but is the best fit to the input 
spectrum. The filled triangles are clearly better matched 
by the dashed line than they are by the solid line. 

However, the higher metallicity synthetic spectrum, 
at [Fe/H] = —0.60, is close to the actual metallicity of 
CD-57:949, as determined from high-resolution analy- 
ses. If the data were more appropriately scaled, as with 
the open triangles, it would better match this spectrum. 
Even in this rough simulation, it is fairly clear that the 
open triangles match the solid line much better than the 
filled triangles match the dashed line. However, from 
the data supplied to autoMOOG, an incorrect, and al- 
ways lower metallicity, synthetic spectrum is chosen as 
the best fit. This problem becomes more severe at higher 
metallicities as the apparent continuum is further sup- 
pressed relative to the true continuum, resulting in larger 
deviations to lower selected metallicities. Thus, this is a 
problem with the input spectra, not with autoMOOG 
itself. With properly continuum-subtracted spectra, au- 
toMOOG will likely produce results with similar scatter 
seen at low metallicity. 

Efforts have been made to include renormalization 
within autoMOOG to attempt to deal with this prob- 
lem. While renormalization does show some promise, 
its results are currently erratic. Future versions of au- 
toMOOG will likely include renormalization, but it will 
take considerable modifications to the code before this 
becomes a stable feature. 

The line index results (Figure [5]d) exhibit more scat- 
ter than autoMOOG at low metallicity, but similar be- 
havior at high metallicity. For the line index approach 
at high metallicity, this is caused in part by the same 
continuum issues mentioned above but also by satura- 
tion of the KP index. This can occasionally be mini- 
mized by substitution of other metallicity determination 
methods designed specifically to handle higher metal- 
licity star s, such as the auto-corr e lation function sug- 
ge sted by iRatnatunga fc Freeman! (|1989( ) and explored 
bv lBeers et al.l ( 19991) . This accounts for the smaller de- 
viations from the high-resolution results of the BM stars, 
but not the HERES stars, seen in the line index compar- 
ison plot. Also note that the stars with [Fc/H] < —1.0 
that show the largest deviations in Figure [3^, arc simply 
not present in Figure [3Jd, as the line index method failed 
to determine any metallicity for these stars. 

As such, metallicities determined with autoMOOG are 
generally as good as those determined using existing 
techniques, and somewhat better at the lowest metal- 
licities ([Fe/H] < —3.0). Since this method was designed 
for use with VMP stars, it is then better than existing 
techniques at accurately determining metallicities. How- 
ever, due to the difficulties at higher metallicity, it should 
not be used blindly with samples that may contain stars 
with metallicity greater than about [Fc/H] = —2.0. 

Much better agreement is seen when comparing the au- 
toMOOG fits on carbon, [C/Fe], to the high-resolution 
abundances (Figure 0. The only stars plotted are those 
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Fig. 5. — Comparisons of the carbon abundance, [C/Fc], deter- 
mined from autoMOOG to that determined from high-resolution 
analyses, plotted vs. the high-resolution (a) [C/Fe] and (b) [Fe/H]. 
The points are the same as in Figure [3] Upper limits are not plot- 
ted, only determined values. 



for which fits (as opposed to upper limits) have been 
obtained and for which high-resolution abundances are 
available. When plotted as a function of [C/Fe] (Fig- 
ure [5^), it is seen that there is very little scatter at 
lower [C/Fe] and only slightly larger scatter as the carbon 
abundance is increased. At higher carbon abundances, 
the added scatter is likely due to the presence of other 
strong carbon features in the vicinity of the G band, or 
the decreased sensitivity of the G band as it begins to 
saturate. 

Examination of the carbon fits as a function of metal- 
licity (FigurcEb) indicates little scatter, but a small slope 
seems to be present. This slope is not real, but is caused 
by a jump at [Fc/H] ~ —2.5, with a slight underestima- 
tion of [C/Fe] at lower metallicities, and a slight over- 
estimation at higher metallicities. The underestimation 
at low metallicity is likely caused by the use of differ- 
ent carbon features at high and low resolution, which 
are affected differently by model atmosphere parameters 
(i.e., one-dimensional versus three-dimensional and LTE 
versus NLTE). This then should extend to higher metal- 
licities as well, for a global underestimation of around 
0.25 dex. 

The slight overestimation of [C /Fe] at high metallicity 
is likely caused by the same difficulties that have already 
been seen at high metallicity, although the effect is much 
weaker here. The reason for this weaker effect is that 
the carbon fits are actually a fit of [C/M], where M is 
the scaled-solar model abundance. To ultimately derive 
[C/Fe], this model-dependent carbon value needs to be 
offset by the deviation of iron from the model. This 
offset is at most half of the distance between models in 
our grid (i.e., A [Fe/H] < 0.25 dex), and thus should not 
significantly affect the final results, but causes this slight 
relative overestimation. 

In addition, the temperatures used were determined 
from photometric colors. At times, the parameters 
adopted for the high-resolution analysis vary from these 
values, occasionally up to a difference of nearly 1000 K, 
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Fig. 6. — Carbon abundances ([C/H]) plotted vs. effective tem- 
perature. Plotted are the abundances for the BM sample (open 
squares) and the HERES sample (filled triangles). All upper limits 
are denoted with an arrow. 

although generally only by a few hundred kelvin. As seen 
in Section 3, such offsets in temperature could lead to off- 
sets in metallicity of ~0.15 dex and ~0.10 dex in [C/Fe]. 
Due to the inconsistent nature of these offsets, this could 
also lead to enhanced scatter in any of the comparisons 
to high resolution. 

It is also possible to use these results to perform an in- 
ternal consistency check on the method. Since a number 
of stars from the BM sample have been observed multiple 
times at medium resolution, the abundances determined 
for each individual spectra can then be compared to see 
if autoMOOG determines the same relative abundances 
for a given star, regardless of the spectra used. For each 
star, an average metallicity and carbon abundance are 
calculated from the available results, and then the resid- 
uals between each individual abundance and the average 
are calculated. 

For both metallicity and carbon abundance, a very 
small scatter around zero of ~ 0.02 dex is observed, well 
within the intrinsic errors, showing the overall stability of 
autoMOOG. The scatter for [Fe/H] at high metallicitics 
may be slightly larger than elsewhere. This is assumed 
to be a result of the exact continuum fit used for each 
spectra, and is thus associated with the general high- 
mctallicity difficulties. Even so, this scatter is still quite 
small, increasing by no more than 0.01 dex. 

5. UPPER LIMITS 

For the purpose of this work, a 67 band is considered de- 
tectable if its deflection from the continuum is at least 5% 
the value of the continuum at that position. Upper lim- 
its for [C/Fe] were determined by varying abundances of 
synthetic spectra and comparing their strengths against 
an essentially "zero" carbon abundance spectrum which 
was used to set the continuum level. A surface in T e ff- 
log g space was then fit to these values. Stars with carbon 
abundances below this surface were reassigned with the 
value of the upper limit at that position, and flagged ap- 
propriately. The approximate position of this cutoff can 
be seen in Figure [HI 

When the 67 band can be properly fit, appropriate 



130 I 

Total Sample — I I 




-5 -4 -3 -2 -1 



[Fe/H] 

Fig. 7. — Metallicity distribution function for the full HERES 
and BM sample (solid line) and those stars for which only an upper 
limit on carbon have been determined (dashed line). 

[C/Fe] values are generally determined. It is, however, 
useful to explore under what physical conditions only an 
upper limit is obtained. The detectability of the 67 band 
relative to temperature and metallicity is shown in Fig- 
ures [6] and El respectively. 

In Figure El the carbon abundance, [C/H], is plotted 
versus effective temperature using the values from auto- 
MOOG. The stars in these samples cover a fairly wide 
range of temperatures, from 4000 K to 7000 K. Both sam- 
ples of stars cover about the same temperature range. 

The G band is a molecular feature, which, for the 
same abundance of carbon, will appear stronger at lower 
temperatures. As temperature is increased, the G band 
weakens until it is no longer detectable. The point where 
this happens is obviously dependent on the level of car- 
bon enhancement, [C/H], and is easily seen in the plot. 
[C/H] is used for this analysis as it is solely a measure of 
the carbon abundance, without any implicit assumptions 
about overall metallicity, which, along with the temper- 
ature, is primarily what the strength of the 67 band de- 
pends on 8 . From the line of detectability, it can be seen 
that for cool stars, the 67 band can essentially always be 
detected. However, for warm stars, especially above 6000 
K, carbon abundances must be fairly high to be detected. 
Indeed, for temperatures above 7000 K, it would be dif- 
ficult to measure a carbon abundance at all, no matter 
how high its true abundance. In part, this is due to the 
apparent upper limit on [C/H]. 

Although there are a handful of stars with slightly 
higher carbon abundances, a general limit to the allowed 
carbon abundances falling at [C/H] ~ —0.5 can be seen, 
which is higher than some ha ve found (—1.0 and —0.7 : 
lRvanll2003t iFrebeT ct al. 20 061, respect i vely), while lowe r 
than others (0.0 : lLucatello et al.ll2005tlAoki et al.ll2007h . 
It is likely that the true location of this upper limit is 

8 The G band also depends weakly on surface gravity. In addi- 
tion, nitrogen and oxygen may lock up some of the carbon in CN or 
CO molecules, which cannot then form CH. However, this is likely 
to produce only a negligible change in the strength of the G band 
at this precision even in cool O- or N-rich stars, where such effects 
would generally be strongest. 
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temperature dependent, with cooler stars exhibiting a 
lower value, and warmer stars a higher value. This may 
be visible in Figure and can further be supported by 
the previously mentioned results, and is presumably due 
to evolutionary effects. 

While the two samples we have examined do not fully 
complement each other, as seen in their combined Mctal- 
licity Distribution Function (MDF; Figure [7|), their re- 
sults can still be used to analyze the detectability of the 
G band with respect to metallicity. Figure [7] shows the 
MDF of the full samples of stars, and the MDF of stars 
with only a measured upper limit. 

This second set of stars overall traces the shape of the 
full MDF, although accounting for slightly larger per- 
centages of stars at lower metallicitics. This is not sur- 
prising, as metal-poor stars have, by definition, low abun- 
dances and weak lines. However, due to the high quality 
of these samples, carbon abundances are able to be de- 
termined for a significant fraction of stars, even at the 
lowest metallicities. 

6. CONCLUSIONS AND FUTURE WORK 

We have developed a new method for the determina- 
tion of [Fc/H] and [C/Fc] for stars using automated spec- 
tral synthesis. This method appears to not be affected by 
the presence of additional C-related feature s in the spec- 
tra of some stars that has been claimed by ICohen et alJ 
(|2005h to limit the accuracy of the line index method 
for CEMP stars, and does not require the massive cali- 
bration required by that method. autoMOOG is signifi- 
cantly quicker and more consistent than spectral synthe- 
sis by hand, which in general could vary from individual 
to individual, and even for one individual over the course 
of time. 

We have explored the behavior of this method for two 
samples of stars with both medium-resolution and high- 
resolution spectroscopy available. Using this method on 
two samples of stars which cover a wide range of param- 
eter space, we are able to see exactly how this method 
behaves. For stars with [Fc/H] < —2.0, both the metal- 
licity, [Fe/H], and carbon abundance, [C/Fe], where de- 
tectable, are able to be successfully recovered. Despite 
the known problems with their metallicities, the car- 
bon abundances for the more metal-rich stars are sim- 
ilarly recovered. As this method was designed and opti- 
mized primarily for use with metal-poor stars, it should 
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not be used indiscriminately with samples containing 
more metal-rich samples, as these stars will contaminate 
metal-poor studies due to the metallicity underestima- 
tion. In future versions of this technique, we will attempt 
to resolve the issues present at high metallicity through 
implementation of various methodologies currently un- 
der development. However, at least for stars with [Fe/H] 
< —2.0, it is more consistently accurate (with errors of 
around 0.2 dex for both [Fe/H] and [C/Fe]) than previ- 
ous abundance methods, especially at the lowest metal- 
licities, where these other methods break down. 

Although we do not show it here, this program is ca- 
pable of fitting any element from any detectable line. 
Thus, it would be possible to, for instance, measure the 
abundance (or an upper limit) of s-proccss elements such 
as barium (4554 A) or strontium (4078 A), expected to 
be enhanced for many CEMP stars. In this way, we 
could roughly separate CEMP stars into the C EMP-s 
and CEMP-no classes (|Beers fc Christlieb! I2005T ) . How- 
ever, detailed studies of these fits have not yet been per- 
formed. In particular, no similar upper limit mechanism 
is in place for these other lines. For the moment, fit- 
ting other features should only be done for very rough 
estimates, and should not be used for any detailed sci- 
ence, but can be used for identification of likely neutron- 
capture enhanced targets for high-resolution follow-up. 
Future versions of autoMOOG will allow fitting for a 
larger variety of lines with more robust detection and 
saturation algorithms, as well as internal fitting of ef- 
fective temperature and surface gravity, to provide less 
dependence on external techniques. 
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